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We demonstrate that a site-dependent driving of a periodic potential allows for the controlled 
manipulation of a quantum particle on length scales of the lattice spacing. Specifically we observe 
for distinct driving frequencies a near depletion of certain sites which is explained by a resonant 
mixing of the involved Floquet-Bloch modes occurring at these frequencies. Our results could be 
exploited as a scheme for a site-selective loading of e.g. ultracold atoms into an optical lattices. 
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Introduction The nonequilibrium physics of driven 
lattices has long been the subject of intensive theoretical 
research [IHi with applications to numerous experimen¬ 
tal setups such as terahertz driven semiconductor het¬ 
erostructures [ 1 , or cold and ultracold atoms loaded 
into driven optical lattices 01. Thereby, it was shown 
that the interplay of a spatially periodic lattice poten¬ 
tial and a driving force leads to a plethora of interesting 
non equilibrium phenomena, a paradigmatic example be¬ 
ing the celebrated ’ratchet effect’ where particles undergo 
directed motion despite the absence of any mean forces 
[iSl. Besides that and triggered particularly by the 
upcoming ultracold atom experiments the inclusion of a 
driving force was used for the renormalization of the tun¬ 
neling rates between adjacent lattice sites llJ or for the 
engineering of so called artificial gauge fields [IJ, [l^ . 
While the main focus has so far been on global driving 
forces that are the same everywhere in space, it has re¬ 
cently been shown how local modulations of the driving 
allow for extensive manipulations of the particles classical 
dynamics leading to phenomena such as a site-dependent 
particle trapping [Ig, [I^ , the spontaneous formation of 
density waves or the emergence of order by the combina¬ 
tion of disorder and driving [isl - l^ . Only very recently, 
this concept of a spatially varying driving has been ap¬ 
plied firstly in the quantum domain , where the trans¬ 
formation of an avoided- to an exact crossing in the Flo- 
quet spectrum as well as the control of asymptotic cur¬ 
rents have been shown and explained as a consequence 
of the local driving. In the present work we demonstrate 
how the site-dependent driving opens the roadway to¬ 
wards a control of the nonequilibrium dynamics of wave 
packets on length scales of the lattice spacing. We pro¬ 
pose how, by exploiting a resonant mixing of Floquet- 
Bloch modes, quantum states which are highly localized 
on specific lattice sites can be prepared in a very con¬ 
trolled manner. The in this way achieved spatially de- 
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pendent occupation of the lattice sites goes clearly be¬ 
yond previously introduced schemes such as ’light shift 
engineering’ [22| where a spatial dependence of the atom 
density was obtained only on length scales of several hun¬ 
dred lattice sites. 

Setup We consider quantum particles in a time- 
dependent periodic potential of laterally oscillating Gaus¬ 
sian barriers in one dimension (Fig. [T](a)). The dynamics 
is governed by the time-dependent Schrodinger equation 
(TDSE) ih-^A>{x,t) = H{x, t)A>{x, t) with the Hamilto¬ 
nian: 


H{x, t) 


2m dx"^ 


CO 

s=—oo 


/ x-s L-ds(t) 

[ -s- 


( 1 ) 


with m, Vo,A and L being the particle mass, potential 
height, barrier width and lattice spacing respectively. 
ds{t) = Acos{ujt -|- ^s) is the driving law of site s with 
frequency w and amplitude A. We will restrict ourselves 
to cases with Ss = Ss+up for some Up €N such that UpL 
is the spatial period of the lattice. 

Formalism and computational scheme We briefly out¬ 
line the employed computational method pfl ]. The tem¬ 
poral periodicity of the Hamiltonian [H{x,t) = H(x,t + 
T) with T = 27r/a;] ensures that every solution of the 
TDSE can be written as A/a{x,t) = t) with 

the real quasi energy (QE) Cq, G [—hoj/2, +Huj/2] and the 
Floquet mode ^a{x,t) obeying ^a{x,t) = ^aix,t + T). 
Due to the additional spatial periodicity [H{x,t) = 
H(x + UpLjt)] each Floquet mode can be written in 
terms of a Floquet-Bloch mode (FBM) as ^a,K{x,t) = 
e"'^'^(t>a,K.{x,t) with (l)a,K{x,t) = 4>a,K{x + npL,t) and 
K G [—7r/(npL),-|-7r/(npL)] being the quasi-momentum. 
For every value of k, the FBMs are eigenstates of the 
time evolution operator over an entire period of the driv¬ 
ing t/«(T + to, io) M- 


C/”(T-bto,io)$a.K(a;,io) = e <o)- (2) 

Hence, by using that the FBMs constitute an orthonor¬ 
mal basis of the Hilbert space spanned by the solutions of 
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with an approximate period of 75T (Fig. [TJb)). Fig.[T](c) 
reveals an additional, seemingly irregular, micro oscilla¬ 
tion of the population ns{mT). In addition, a slow over¬ 
all decay of all three in Fig. (TJc) depicted populations 
is observed, which is a straightforward consequence of 
the diffusion of the initial state. Note that the observed 
oscillations of the population crucially rely on the site- 
dependent driving and disappear when the phase shift 
6 of the central barrier approaches zero (Figs. [Hd) and 
(e)). In the following we will analyze and explain the 
oscillatory behaviour of the site populations naijnT) and 
demonstrate afterwards how it can be exploited for a se¬ 
lective loading of the lattice. 

To understand the origin of the observed phenom- 


FIG. 1. (a) Snapshot of a site-dependently driven lattice 

with lattice site occupations Us- Shaded barriers indicate 
the barriers equidistant equilibrium positions, (b) Time evo¬ 
lution of lattice site populations for an initial Gaussian of 
width cr = 207r in a lattice with 3 barriers per unit cell 
with phases (0,7r, 0). (c), (d) and (e) depict populations 

of the first three lattice sites (encompassed by a black rect¬ 
angle in (b)) for phases of the central barrier in each unit 
cell of TT (c), 0.27r (d) and 0 (e). Remaining parameters are 
L = 10, Vo = 1.0, to = 1.0, A = 1.0, m = 1.0, ft = 1.0 and 
A = 0.5. 
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the TDSE, the stroboscopic time evolution of any initial 
state lE'(x,to) is given as: 

\E'(a:, to + mT) = 

/ -|-■7^/(T^pL) 

dK ^ to). 

-7r/(npL) Q, 

( 3 ) 

where Ca,K{to) are the overlap integrals between the ini¬ 
tial state il>{x,to) and the FBM $a,K(a^,^o)- After nu¬ 
merically diagonalizing the time evolution operator and 
thus obtaining the FBMs as its eigenstates we are able to 
propagate arbitrary initial states by employing Eq. ©■ 
Note that in principle, the FBMs and thus the dynamics 
of an initial state depend on the initial time to- However, 
the findings of this work are very robust with respect to 
variations of to and thus we will only present results for 
the case to = 0 in the following. 

Lattice site occupations for site-dependent driving We 
now demonstrate how a site-dependent driving can be 
used to modulate and control the occupation of spe¬ 
cific lattice sites. We consider a Gaussian initial state 
'h(a;, 0) = (7rtT^)“^/^e“^ that is exposed to a lattice 
with three barriers per unit cell (rip = 3) and Sg = (0, S, 0) 
for s = 0, 1, 2 and with 6 € [0, tt]. The occupation of lat¬ 
tice site s after m periods of the driving is given by: 

psL 

ns{mT) = / \'^{x,mT)'\^ dx (4) 

J(s-1)L 

As a striking feature we observe pronounced and site- 
dependent temporal oscillations of the population Ugit) 


FIG. 2. Setup consisting of three barriers with phases Ss = 
(0, TT, 0) for periodic boundary conditions, (a) Time evolution 
of the lattice site occupation na{mT), (b) and (c) show the 
same quantity but only the 3 (b) or 2 (c) most occupied FBMs 
are taken into account, (d) position representation of the three 
most occupied FBMs. Remaining parameters as in Fig. [T] 

ena let us systematically simplify the setup. As a first 
step, we investigate the population dynamics ns{mT) in 
a setup consisting of only a single unit cell with periodic 
boundary conditions and with a uniform initial state (re¬ 
sembling the conditions at the center of a broad Gaus¬ 
sian initial state in an infinitely extended lattice). Note 
that the inclusion of periodic boundary conditions corre¬ 
sponds to a restriction to k = 0, as every FBM fulfills 
^a,K{x npL,t) = e™p^'^i^a,K{x,t). In fact Fig. [^a) re¬ 
veals that our simplihed lattice is capable of reproducing 
the key features of the population dynamics observed in 
Fig me). Further simplifying the setup, we now perform 
a few mode approximation after sorting the FBMs at 
K = 0 according to their overlap with the uniform initial 
state. Hence, the ’Floquet ground state’ (FGS) corre¬ 
sponding to the largest overlap is labeled as $o, the mode 
with the second largest overlap as $1 etc. Evidently, for 
the three mode approximation (i.e. we set Ca = 0 for 
a > 2 in Eq. the pronounced oscillation of the lattice 
site population is reproduced (cf Eig. HKb)), while it is 
not for the case of only two modes (cf Fig. HKc))- This 
suggests that the dynamics of the lattice site occupation 
is predominantly determined by the FBMs $0 and <I> 2 , 
as $1 does not contribute to the pronounced oscillations 
(cf Fig. me) and see discussion below). Let us therefore 
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assume that only these two modes are initially occupied, 
i.e. 0) = C'o$o(2;)0) -\- C 2 ^ 2 (.x,Q). According to 

Eq. [3] we obtain: 

|«'(x,mT)p « |Co$o(a;,0)p + \C2<^2{x,<^)\^+ 
2Re(C'oC'2*$o(a:, 0)$;(a:, 

(5) 

The interference term oscillates with a period of Tqsc/T = 
Sw/(e 2 — Co) ~ 77, which fits the period of the population 
oscillations in the infinitely extended lattice (Fig. [Ijb)) 
quite well. A remaining question is why a few mode ap¬ 
proximation which includes $2 leads to pronounced os¬ 
cillations of ns(mT), but not the corresponding approx¬ 
imation based solely on $0 and $ 1 . The reason for this 
is that the interference term in Eq. [5] can only lead to 
a pronounced oscillation of UsimT) if at least one of the 
involved FBMs is strongly localized on one of the lattice 
sites. As Fig. [5] (d) reveals, this is indeed the case for 
$ 2 , but not for $ 1 . Hence, even though the FBM $1 is 
more populated than $ 2 , it contributes much less to the 
distinct time evolution of the observable ns{mT). 
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FIG. 3. (a) Population of the s = 0 lattice site for a setup 
containing three lattice sites with phases (0, tt, 0) and periodic 
boundary conditions {k = 0) for different driving frequencies 
u. (b) Extract of nniax(ta) (black) together with the overlap 
of the FGS with a uniform state 'Fu (blue). Vertical lines are 
by Eq. dS]) predicted resonances for n = 1 (solid) and n = 2 
(dashed). (c) Maximal population of any of the 3 lattice sites 
nraax((jj) which is reached within the first 400 driving periods. 
Inset of (c) shows the QE spectrum (small arrow indicates the 
FGS). Remaining parameters as in Fig. [T] 

Controlling the lattice site populations We now 
demonstrate the controllability of the time-dependent 
populations of the individual lattice sites via the driv¬ 
ing frequency w. First, we keep the periodic boundary 
conditions, i.e. k = 0, and discuss the impact of nonzero 
K afterwards. In the low frequency regime (w < 3.5), the 
population of site s = 0, no(mT), exhibits oscillations 


whose period depends sensitively on oj thereby ranging 
from Tosc 5T to Tqsc > 400T (Fig. [3Ka)). Contrarily, 
for u! > 3.5, no{mT) becomes approximately constant 
with nQ(rnT) 1/3 except at very narrow frequencies 
intervals (e.g. at w ~ 5) where n^lrnT) again strongly 
deviates from 1/3. Since the population dynamics of the 
sites s = I, 2 is qualitatively similar we refrain from show¬ 
ing it. Instead, we show the maximal occupation rimax of 
any of the three lattice sites that is reached during the 
first 400 driving periods as a function of w (Fig. [3Kc)). 
The regime w 3.5 exhibits an irregular behaviour re¬ 
flecting the above observation of a high sensitivity of the 
lattice site populations on w. Besides narrow peaks, we 
observe for w > 3.5 that nniax(w) is approximately 1/3 
corresponding to an equal population of the three sites 
at all times. The overall convergence for w —00 of 
ns{t) towards 1/3 is not surprising, as for fast driving the 
Hamiltonian can typically be approximated by its time 
average, which is site-independent in our case. Thus, it 
remains an intriguing question why the occupation of the 
different lattice sites inherits such a strong imbalance at 
distinct driving frequencies even in the regime of large w. 
To gain insight, let us recall that the oscillations in the 
population dynamics ns(rnT) can be attributed mostly 
to the interference of the FGS $0 with one of the higher 
modes. Since $0 is the dominantly occupied FBM, a 
strong desymmetrization of the lattice site populations 
as indicated by a large value of rimax can be expected 
to correspond to a large desymmetrization of $ 0 . Mea¬ 
suring the desymmetrization of $0 by its overlap with a 
uniform state = 1/ \/npL, we find that the peaks in 
’T-max(w) are indeed accompanied by resonant dips of the 
overlap of $0 with (Fig. [3Kb)). 

The remaining open question concerns the physical mech¬ 
anism underlying the observed desymmetrization of $0 
at specific frequencies. Here, the frequency-dependent 
QE spectrum reveals further insight (Fig. jSKc) inset). 
While the spectrum appears irregular at small oj < 2 
we see clearly the emergence of QE bands at larger w. 
Moreover, these bands are subject to multiple exact- and 
avoided crossings. Of particular interest to our analy¬ 
sis is the FGS $0 (marked by the arrow in Fig. |3Kc) in¬ 
set). A closer inspection reveals that each of the peaks in 
?T'max(w) and equally the dips in | (d’ol'F^) ^ fall together 
with avoided crossings of $0 with states belonging to 
higher bands. Hence, the FGS mixes strongly with an¬ 
other state at these frequencies which leads to the strong 
desymmetrization and thus to a resonant enhancement 
of rii„ax(w). By exploiting this insight we can employ a 
remarkably simple approximation in order to deduce the 
resonance positions in frequency space. As FBMs with 
average kinetic energies ^ Vq represent an almost free 
propagation through the lattice we approximate Vq ~ 0 
and calculate their QEs easily as Ca = 2'K^a^/m{npLY 
(which then have to be folded back to the first Brillouin 
zone). In this regime of fast driving we may further as¬ 
sume eo ~ 0 compared to the extension of the Brillouin 
zone huj. Gonsequently, the frequencies Wres correspond- 
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ing to resonances in nniax(w) caused by avoided crossings 
of the FGS and the higher bands, can be approximated 
as the zeros of the free particle quasi energy: 


a;res(Q!,n) 


2'K^h 

m{npL)‘^ n ’ 


for n = 1, 2, 3,... 


( 6 ) 


Here, n denotes the number of times the free particle 
energy has to be folded back into the first Brillouin zone. 
In Fig. |3jb) we show the locations of the resonances 
predicted by Eq. m for n = 1,2 and up to a = 20, 
corresponding to the number of FBMs which are taken 
into account numerically and which are shown in the 
spectrum (cf. inset Fig. [3] (c)), confirming a good 
agreement. 



FIG. 4. (a) Time evolution of the lattice site occupation 

ns(rnT) for the resonant case w = 2.74. (b) ns{mT) for 

s = 0,1,2 (encompassed by a black rectangle in (a)), (c) 

same as in (b) but restricted to k = 0. The s = 0 and s = 1 
curves in (b) and (c) are barely distinguishable. Remaining 
parameters as in Fig. [T] 


Resonant population transfer at nonzero quasi¬ 
momentum Finally, we demonstrate how the above an¬ 
alyzed resonant population imbalance is altered by the 
inclusion of higher quasi-momenta. We explore the time 
evolution of the same Gaussian wave packet as above 
(Fig. [IJb)) but now with a driving frequency uj = 2.74 
where nniax(w) features one of its strongest peaks. After 
a transient of t ~ 50T, we observe a pronounced imbal¬ 
ance of the lattice site populations reaching its maximum 
at t/T m 250 (FigH] (a) and (b)). Interestingly and in 
clear contrast to the k = 0 case (Fig.HKc)), we observe no 
indication of a pronounced oscillation of population be¬ 
tween the sites. Apparently, the inclusion of higher quasi¬ 
momenta leads to a decay of the oscillation amplitudes. 
Thus, while we still observe an overall diffusive decay at 
large times, the imbalance between the lattice sites is ef¬ 
fectively frozen at intermediate times (500 < t/T < 800). 
Note that in the previously studied case of a; = 1 (Fig [1] 
(c)) we could not see this effect as the oscillations there 
occur on much faster time scales compared to the de¬ 
cay of the oscillation amplitude. To gain further insight, 
let us make some simplifying assumptions. According 


to our previous analysis, the dynamics at k = 0 can be 
reasonably well described by including only the FGS $0 
and one of the higher modes, say $e. Gonsequently, we 
may obtain some intuition on the main features of the 
K ^ 0 dynamics by considering the two corresponding 
QE bands: $o,k and ^e,K, so that the time evolution of 
4'(a;,0) reads (cf Eq. [3]): 




-|-■7r/(r^pL) 


—■7r/(npL) 




a 






(7) 

Numerically, we find that the periodic part of the FGS 
4 )o,k and the associated QE depend only very weakly 
on the quasi momentum k, i.e. ~ and 

eo.K ~ Co, which is in contrast to the excited state 
and ee,K respectively. Thus, upon evaluating 
we encounter a term proportional to |4)o(a;,0)p plus 
terms including the k integral over .K{a:,0). 

The latter ones, will quite generically inherit a diffusive 
decay due to the k— dependence of Ce^K and thus tend 
to zero for large mT. Gonsequently we are left with 
|4>(a:,TOr)p oc |<l>o(a;, 0)p and the occupation of the lat¬ 
tice sites becomes indeed time-independent and simply 
follows the lattice site occupations as given by the FGS 
4)0(a;, 0) which features a strong desymmetrization be¬ 
tween the lattice sites. Note that this simplified line of 
arguments neglects -among others- the diffusion caused 
by the weak K-dependence of 4 )o,k or by driving induced 
couplings of <i)o,K to higher bands which lead to the over¬ 
all decay even at late times as seen in FigH] (b). 

Conclusion Our study of the temporal evolution of 
the lattice site occupations for a quantum particle in a 
site-dependently driven lattice demonstrates that we can 
achieve an almost complete depletion of certain lattice 
sites by exploiting narrow avoided crossings of the Flo- 
quet ground state with higher Floquet bands. This rep¬ 
resents a mechanism for a controlled site-selective load¬ 
ing of a lattice potential, which might open the road¬ 
way towards the preparation of initial states which are 
otherwise hard to obtain. The latter are certainly of 
interest for e.g. the investigation of the highly corre¬ 
lated many particle systems in optical lattices and their 
follow-up nonequilibrium dynamics. In ultracold atomic 
physics Feshbach resonance management could be used 
to appropriately tune the interaction strength among the 
atoms [24| after the preparation of the corresponding ini¬ 
tial state. Experimentally, our results are of relevance 
whenever the translational symmetry over at least two 
neighbouring lattice sites is broken. In particular, for 
cold atom experiments performed in shaken optical lat¬ 
tices, there are several ways how this could be imple¬ 
mented, among them Fourier synthesized lattices where 
higher order Raman transitions allow for a substructure 
below half the lasers wavelength [ 2 ^, distorted 2D lat¬ 
tices [ 2 ^ or so called ’painted potentials’ where full con¬ 
trol over the potential landscape is achieved 
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